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ABSTRACT 


In this thesis, a new model is suggested for 
modelling signal data from a deterministic or stochastic 

process. The signal data is modelled as sum of several 

complex Frequency Modulated Signal. Autoregressive model 

based on Burg algorithm is used for estimation of model 

coefficients. These coefficients are used to find the power 

spectral density. From this the FM subsignals y Carrier and 

Modulating Frequencies are estimated. Use of Discrete 

Fourier Transform is made of in estimating the subsignals 

Modulation Index. Separation of individual subsignal before 

estimation of Modulation Index is also discussed. These 

estimated parameters are then used in estimation of 

subsignals amplitude and Phase. 

The model is first fitted on a synthesised data. 
The model is then fitted on speech signal with varying 
degree of success. 
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CHAPTER- 1 


INTRODUCTION 

Parametric model I ling of signal is an important step in many 
signal processing applications- Many practical signals, like image 
and speech signals , are highly redundant and it may be desirable 
to represent the signal with a fewer number of samples for economy 
of storage and /or transmission bandwidth limitations. Knowing the 
parameters of a signal the resulting parameters can be used as a 
feature vector for pattern recognition , speech recognition and 
alanguage identification etc. 

A large number of signals show predominant Frequency 
Modulation in the time domain. As such they can be assumed as 
consisting of a combination of several complex Frequency Nodulated 
< FM) singnals, eg speech signals etc. In this thesis an attempt 
IS made to decompose a signal as a combination of various 
frequency modulated signals. 

1.1 NON- STATIONARY SIGNALS 

The term stochastic process or random process is used to 
describe the time evolution of a statistical phenomenon according 
to probabilistic laws. One particular realization of a discrete 
time stochastic process is called a time series. For example , the 

sequence x (n), x (n-1>, x(n-N+1) represents a time series 

that consists of present observation x(n) and N-1 past 



observation of the process made at time n— 1, n-M+l .We say 

that a stochastic process is strictly stationary if its 
statistical porperties are invariant to a shift of time , le the 
Joint probability density function of these observations made at 
times n, n-1, n-h+1 must remain the same no matter what 

values we assign to n for fixed M. 

In practice, it is not possible to determine the joint 
probability density function for an arbitrary set of observation 
made on a stochastic process. As the assumtion of strict 
stationarity is too restrictive, a notioin of wide sense 
sationarity ( W.S.S) , is used demanding only the invariance of 
the first two moments with time. A discreate time stochastic 
process is said to be wide sense stationary or weakly stationary 
if its mean is constant and its auto — correlation depends only on 
the lag le. 

E Cx (n)3 = (j and E Cx (n) x*( n + k)D = r( k ) ; k = O, ±1,±2.. 
Where E = Expectation operator 

» = Complex conjugate operator 

Any signal which is not wide sense stationary is considered 
as a non-stationary signal C 1 D . 

1. a. PARAMETER ESTIMATIOM OF NON-STATIONARY SIGNALS 


Parametric 


method for analysis and synthesis 


of 



stationary signals is a well developed field .Also well developed 
are Nonpararoet ri c methods of nonstationary spectrum estimation. A 
lot of work has also been done in modelling non-stationary signals 
within varying degree of success .Adaptive filters have been 
extensively used for modelling non-stationarity » processes. 
However, adaptive filters are useful, only if the measure of 
non—stationarity , which may be the rate of variation of 
autocovariance function with respect to time is less than the 
speed of adaptation C S 3 . If this condition is satisfied , then 
the error in estimating the weight vector has two components .The 
first one, called the weight vector noise, is proportional to the 
speed of adaptation, and the second component called the weight 
vector lag is inversely proportional to the speed of adaptation. 
Hence, one optimum operating point is chosen to keep the total 
error minimum. The conflicting requirements of the components of 
error , limits its use in the case of highly non-stationary 
process. 

Many deterministic and stochastic discrete time processes 
encountered in practice are well approximated by a Rational 
Transfer Function (RTF) model. The RTF model converts the signal 
representation problem into an equivalent parameter estimation 
problem Autoregressive (AR) 5 Moving Average (MA) and 

Autoregressive Moving Average <ARMA) models have been successfully 
used for stationary signals (1^ ) . Application of these methods 

to non-stationary signal requires initial processing of signal so 
as to have a stationary signal . Further , peak of the PSD though 
occurs at the signal frequencies, but its magnitude is not 



proportinal tothe power of the sinusoid as is the case for a 
Fourier based spectral estimators. The thesis lay out is as 
described in the following paragraphs. 

Chapter 2 deals with the method for frequency determination 
for the non— stationary model under consideration using the PSD 
approach. It also discusses the use of Discrete Fourier Transform, 
for estimation of parameters influencing the power of the signal. 
A method for resolving overlapped signals in a low point DFT is 
also discussed. 

Chapter 3 deals with the application of the theory developed 
to synthetic data in cases where the carrier frequencies are wide 
apart and where they are close. 

Chapter 4 deals with the application of the model to speech 
signals with varying degree of sjifccess. 



CHAPTER - 2 


2.1 FREQUENCY MODULATION 

Modulation may be generally defined as the process 
whereby a message is transformed from its original form, into a 
signal that is more suitable for transmission and processing to 
meet needs imposed by specific circumstances. Mathematically, 
modulation may be described as the process of mapping from a 
message space to a signal space. 

The realisation of the sensitivity of amplitude 
modulation to noise, interference and distortion lead to the 
conception of Frequency Modulation (FM) . Frequency modulation is 
a process in which an information bearing of signal , g{t) is used 
to vary the instantaneous frequency©^ another signal c{t>. The 
signal g(t) is referred to as the modulating signal and the signal 
c{t> as the carrier signal. FM belongs to the class of non-linear 
modulation unlike AM which belongs to class of linear modulation 
C6D. 


Consider frequency modulation of a complex exponential 
carrier by a complex exponential modulating signal, so that 
c(t) » A^e-’^^c^ * ^c^ and g(t) = A^e-^ 



The frequency modulated signal xCt) is given by ? 

= c<t> 


yv j {ft t + <^c) j kX A 
x{t)=A^e c ^ -e a 


ilTl t 

^ j Cftct <^c) ^ e m 


(H-1) 


where ft = S"” ~ ~ modulation index . More precisely 

m m 


x(t) can written as 

^ ft cos ftt j<ftt+«* + /55in ft t) 
x(t)=A^e me c ^c m 


(S-S> 


This equation clearly shows that the modulating signal not 
only affects the instantaneous frequency of the carrier signal but 
also its amplitude. That is to say Eqn <2. 2) represents an Am-FM 
combine case. 


For only FM consideration we modify the Eqn (2.2) in the 
following manner 


g(t) = A cos ft t 
^ mm 


Then 


K<t) = 

= A eJ **^*P sinOJ) 


(2-3 


A discrete time version of Eqn (2.3) will be given by 


xcna = A^ e-’'“c" 


n * 0,1, . . . , N-1 (2-4 

Where <ji - OT is the discrete frequency with T as the sampling 


interval 



2.2 SIGNAL MODELLING 


We will consider here a signal model comprising of a 
superposition of several frequency modulated subsignals. A complex 
signal x(t) consist— ing of P complex exponential carriers with 
sinusoidal modulation can be expressed as 


x(t) =1;; A e^ ^^1^ "^1 ^1 

1=1 ^ 


< 2 - 2 . 1 ) 


where = amplitude of i*" exponential carrier 


^1 = phase of 1 exponential carrier, which is a random 
variable uniformly distributed ( 0 ,az). 

— modulation index of 1^^ exponential 
Oi = frequency of exponential carrier 

Ni = frequency of i^^ baseband singal 


Satisfying the sampling theorem, an appropriate sampling 

frequency with time period T is chosen for the uniform sampling of 

x(t). The discrete— time samples of the singal are of the form 

«CnD . E A aJ '“i " * "i"* (e.Z. 2 ) 

i =1 ^ 


where w^ 


i discrete angular carrier frequency 
^ h 

i discrete angular modulating frequency 

i_ «. th 

phase of 1 exponential which is a r.v uniformly 


distributed (o , Sin);- 



The purpose of the thesis is to perform anlysis of a given 
complesi signal so that it can be expressed as a sum of several FM 
subsignals with complex exponential carriers and sinusoidal 
baseband singals. In otherwards» the parametric modelling of the 
signal xCnl will be achieved by determining the parameters of 
equation (2>2 m 2)> The steps involved will be s 

1 . Estimation of the carrier and modulating 

frequencies . 

E . Estimation of the modulation index of each 

sub-signal . 

3 . Estimation of the amplitude and phase of 

each sub-signal . 


2.3 MODIFIED AUTOCORRELATION FUNCTION 


For the signal xCnl which is a sum of complex FM signals, the 
autocorrelation function is given by 

r Cn.kl = E C X* CnI xCn+k3 3 {2.3.1) 

X 


where k is the lag 

from Eqn (2.2.2) we get 
* P -j 

x Cn3 = E A e 

isa-l ^ 


n ^ + rt. ain n) 

1 1 i 1 


(2.3.2) 


P 

xCn+k3 E ^ 

i=1 


j to (n+k)-hi -ifi sin(v . (n+k) ) ) 
e 1 11 i 


(2.3.3) 

from Eqn (2.3.1), (2.3.2), and (2.3.3) we have 



(2.3.4) 


■ Cn,k3 = E [ E A. e 
^ 1=1 ^ 


+ ft^sin v^n) 


E rT| 
1 = 1 ^ 


^ (6ij(n+k) ■*■ sin 




= 0 


o 


E [ e-l ‘^l^l’] , = E C13 

f'l 

In view of Eqns (2-3.5) and (2.3.6) , Eqn (2.3.4) becomes 


(2.3.5) 

(2.3.6) 


P o - I \ f (En+k). 

r Cn,k3 = E A^. =*" ‘*'i'2'> ‘"i ~r ’ 

“ T=1 * 

(2.3.7) 

Because of its double dependence on position n and lag k, the 
autocorrelation function, r Cn,kD can not be used in 

determination of the frequencies (carrier and modulating) of the 
signal. Also, since discrete values of xCnH are available for only 
one observation, ensemble average can not be taken. 


To overcome these shortcomings of the time-dependent 
autocorrelation function, we take the product function 
P CkH, defined as 


£kl = £ X Cnl K Cn+k3 ] ^ 1 =-™ 
* 2 


P^Zk3 = [ E A. e 


k - O, ±2, ±4, ±6,' 

-j {6i> n + + ft. sin V n) 


1=1 


(2.3.8) 



I fl ^jto,<n+k) /JjS.n (n + k)] g , 

1 = 1 ^ 


E E A A, eJ ‘"l" "“i" e-l ‘<^1 "<‘i’ 

1=1 1=1 ^ ^ 


j (^,sin (i»,(n+k)) — ft sin (w.+n)^ 
» 1 1 1 1 


<S.3.10) 


To have further insight of Eqn (2.3.10), we deal with the 
specific case of say p=2. then 

i^lk 

P^CkT = a! 2/5iSin — — ) 

#C I 

i>„k 

^ CWgk + — 2 ^ 


^ ^ \ 4. ^ sin 

+A^Age-’ S a. 1 'i 


(i>1-g) + /^gSin (Wg 


+ A ^ A 2 e - 


- 4>. 


.) + 


/ 9 ^sin 


/ 5 gsin ( t >2 


Let us a alyse first term of Eqn (E.3.11) 


A ^ 






V ,k 
sin 1 


(E.3.11) 


(E.3.11a) 


Using the properties of bessel functions we know CE3. 





(2.3.12) 


^jxsin ^ ^ 


l=-eo 


( x } 


Uhere coefficients Jj^ are Bessel functions o^ the first 
kind, of order 1 and argument x« Hence Eqn (2.3.11a> can be 
rewritten as 


2 -J (».k) ” , 


m 

E 




(2.3.13) 




E <2^^> ^ 2 \ k=0, ± 2, ± 3,., 


(2.3.13 




Similarly the other terms of Eqn (2.3.11) can be expessed in 
terms of the Bessel function, giving the expression of P Ck3 i.e. 

X 

Eqn (2.3.11) 


a> ^ 1^ 1 

P^Ck3 = A^E J (2/? ) e-» ‘^1*^ "■ "“2“ > 

1 HI SSmiWWI. ^ 


1 1=-Tao 


06 


4 •*12 <^2> 

12= -00 


lo 1*0 •< 

.J ("gk -S-g§— ) 


♦ fl e-'‘*2'^l’ £ 

l1=-oo 


oo 


-'ll •'»2 '^h’ 

6> + 6> 1.1>. + 1^ 


lE=-t» 12=-tB * * • At «; 2 

+ (— J-i— 2_2j^ j 

The fallowing observations can be made from Eqn (2.3.14) 


(2.3.14) 



( 1 ) 


The spectrum of P Ck3 will contain the carrier components 
occuring at twice the original values as k takes only even 
values and an infinite set of true modulating frequencies 
located symmetrically on either side of the carriers at 
frequency separation ofw , 32 p , 3^ .... 

(2) The amplitude of the various component varies 
with modulation index. 

(3) Spectrum will also contain cross carriers which are the mean 
of carriers and cross modulating frequencies . 

The magnitude of Bessel functions J^(^) become very 

small for jnj > \ft\ w C63. For all practical purposes jnj = jf5+1|, 

is the largest n which will have any appreciable value for 

J_(/9).C63 . Therefore in the summation of Eqn (2.3. 14) for 
n 

practical purposes only those terms for (n) extending from zero to 
little greater than ft need be included. 

2.4 ESTIMATION OF CARRIER AND MODULATING FREQUENCYES 
2.4.1 POWER SPECTRAL DENSITY TECHNIQUES 

The power Spectial Density (PSD) of a complex WSS random 
process xCn3 is defined as 

E r exp (-jBifk); s! f £ (2.4. 

k=-o& ^ ^ 

The frequency (f) may either he thought of as the fraction of 
the sampling frequency used in obtaining the data samples from a 
continuous process or as the number of cycles/sample. The PSD 



function describes the distribution of power with frequency of the 
process under consideration. Since the PSD is a function of an 
infinite number of ACF values, the attempt of estimating the PSD 
based on a finite data set is errorneous. At best we can only 
estimate a limited number of ACF values, and therefore, parametric 
modelling is employed to extrapolate the ACF sequence for all 
lags. 


This allows us to replace spectral estimation problem by a 
parameter estimation problem. It is critical, however that the 
model be an accurate representation of the underlying process, at 
least as far as the PSD is concerned. 

Thus the spectral estimation, in the context of modelling 
becomes a three step procedure. The first step is to select a 
model. The second step is to estimate the parameters of the 
assumed model using the available data samples. The third step is 
to obtain the spectral estimate by substituting the estimated 
model parameters into the theoretical PSD implied by the model 
USD. 

Many deterministic and stochastic discrete time processes 
encountered in practise are well approximated by a rational 
transfer function model. In this model input driving sequence uCnD 

and the output sequence xCnD that is to model the data are related 
by the linear difference equation 



(2.4.S.S) 


P <1 

KCn3 « - r a Lkl x Cn-kl + £ bCk3 u Cn-k3 
k=1 k=0 


The most general linear model is termed as an ARMA model. The 
trans-fer -Function H(z) between the input uCnT and the output xCnU 
for the ARMA process is the rational function 


£ bCk3 z 


H(z) 


Blz2 

A<z) 


k=0 


£ aCkH z 


k=Q. 

Often the driving process is 
of zero mean and variance 
becomes 


assumed to 
The PSD of 


be a white noise 
the ARMA output 


sequence 

process 


'^ARMA 


(f ) 



Bif) ,2 

A(fT » 


whre A{f) = A<z> 
B<f) = B(z) 
If all the aCkD 
parameters, then the 
and 




j2 nf 
j2 nf 


coefficients, except aC0Il=1, vanish 
process is strictly a MA process of 


for ARMA 
order q 


^MA == "" |B(f)p 

If all the bCkO coefficients except bC03=1, vanish for ARMA 

parameters the process is strictly an AR process of order p and 

2 

^AR = 2" 

jA(f)j^ 

The AR model is appropriate for spectra with sharp peaks but 


not deep valleys. The MA model can be used for spectra with deep 
valleys but not sharp peaks. The complex FM signal model spectra 


4/i 



will yield many pairs of carriers and carrier plus modulating 
frequencies. The AR model therefore appears to be the most 
appropriate for this model. 

The selection of model order in AR spectral estimation is a 
critical one. Too low an order results in a smoohed estimate, 
while too large an order causes suprious peaks and general 
statistical instablity. 

For the compleK FM signal model the product function (PCkTj 
k— 0, ± S, ± 4, ... ) are fitted onto an AR model based on Burg 
algoritnm [1311. The AR parameters so obtained are used to compute 
the PSD of the product function from Eqn (2.3.14). The PSD should 
consist of several sets of peaks, each set consisting of a single 
carrier freqency which is twice the carrier frequency followed by 
true modulating frequencies symmetrically , distributed about the 
carrier with diminishing magnitudes. Also present will be peaks 
representing cross term carrier which are sum of the true carrier 
frequencies followed by the cross modulating frequencies. If the 
data is real then the PSD will be symmetrical about frequency f=0 
and accordingly the AR model order chosen should also be twice the 
number of frequency pairs expected. 

2.5 ESTIMATION OF MODULATION INDEX 

From Eqn (2.3.14) it is clear that the modulation index ft, 
for a given set of carrier and modulating frequencies and lag will 
affect the magnitude of the PSD at those frequencies. The 
magnitude however, is not proportional to the power of the 



sinusoid as in the case for a fourier based spectral estimator but 
to the square of power. For a single complex sinusoid in noise 
with SNR & and model order p, in the limiting case for ^np >> 1 

and p >> 1 

Pan (2.5.1) 

Hr? max w 

The minimum variance spectral estimate though having a 

resolution poorer than AR spectral estimate has magnitude directly 

proportional to power of the sinusoid for a given model order 

2 

le. P^ygg (f) « A (2.5,2) CKAYD 

Both Eqn (2.5.1) and (2.5.2) clearly depicts that the 
magnitude of the spectral estimate is dependent on the model 
order. We utilise the fourier transform approach to evaluate the 
modoulation index. 

Given a complex valued data sequence xCnT of length N, the 
Discrete Fourier Transform is given by 1153 

f^-1 _ arkn 

XCK3 = £ xCn3 e ■* N~ ? O £ 

n=0 

Substituting value of xCnl from Eqn 

( 2 . 2 . 2 ) 

p 

XCK3 = E { r A. e^^'^'i" SI" ^ atkn 

^=0 ^=1 ^ lie ----- 

Let eJ^^i 

From eq (2.3.12), 


K £N-1 
(2.5.1 


(2.5.2) 


we can write 



^j/5sin«n^“ J ,(,?) 

l=-t» ^ 

Hence eq (2.5.2) can be written as 


N-1 p 


ce> 


XCK3 = E < E fiji = ■’ 

n«0 1=1 li 


N 


(2.5.3) 


Rearranging the suiiuiiation signs 

p CO N — 1 

= E A . r J ,.{/?.) E e 

T ^ Cl T li 1 „ 

1*1 li=-<oo n=0 


j {6ii — an k 


N 


+ Iw )n 
1 


(2.5.4) 


M-1 , a nk . , , ^ J )M 

But E .•■‘"i - "S — ♦ i — 

n=0 l-e-* ‘"r . ) 

N 1 


Hence 


XCKD = 


p CO . j(w.— + Iw ) N 

r A . r j,,(^ ) i rj ( 2 . 5 , 5 ) 

Cl ll ‘ 1 L j (w. - 2nk J 

1=1 li-na> 1-e 1 -r^ + Iw 

N 1 


or 

j XCK3j 


P « .l-uJ'"!" * ‘"i*” 

= E A. E J, (/? ) -----yi i 

11 1 ' . j(w.— aik . - 

1=1 li=-<» 1-e 1 — n — + Iv . 

N 1 


(2.5.6) 


Magnitude of the DFT ie |X{IKI}| will have peak values at those 

where for li = O, ± 1, ± E.... 

w — — + Iv . =0 

1 N 1 

or = (w. + li^^) ; 1=0, ± 1, ± 2,.... 


Expanding Eqn (2*5.6) as 
00 


0 K ^ N-1 
1 


CO 


I X CkDj= A E J {ft i r + A E r_ , + 

1 =-® 1 1 11 2 l=-tB 5 zl 2 2,1 


(2.5.7) 


00 


A 


JE J, (/5 ) r - 
Pi=^ 1 p P'l 


(2.5.8) 



where r 


1—6 1 


aik 


N 


+ ly.)N 
1 


1,1 


. j (w . 
i— e 1 


'm 


t—rs 

+ 


At the 1st highest peak corresponding to =* 

1 X '=‘<31= -i- Ag •‘a.o*-” 

Assuming that the values 
sufficiently small, their contribution towards { 
neglected. Hence Eqn. <2.5.9) can be rewritten as 

I X 1= A, (,S^) ^ 

Similarly at 1 = + 1 and -1 
I X Ck23|= A^ </?^) r^ ,, and 


O, 

+ A J, <^ ) K„ o 
p Ip ' p P,0 

(2.5.9) 


, . J- </9 ) are 
Ip p 

X CkD] can be 


C2.5.10) 


(2.5.11) 


j X Ck23|=|A^ r^ _^1 

= (/5^) r^ (2.5.12) 

Taking the ratio of Eqn (2.5.10) and (2.5.11) we get 


A, J^(/5,) 


X_CK1D 
X Ck2D‘ 


or 


>lo</^1> 




1,1 X Ck U 
1,0 X Ck23 


Similarly we can get 




and 

1 p 


Once the ratio of J^(/9') /J ^(/9) is known we can get 
of /9, the modulation index from the tables having 


Jq(k)/J^(x) for different values of (x). 


(2.5.13) 


the value 
values of 



2.6 ESTIMATION OF AMPLITUDE AND PHASE 


Once the true carrier and modulating frequencies and 
corresponding modulation index are determined by the method 
discussed above, the remaining probelm of estimating the 
amplitudes and phases turns out to be a linear estimation problem. 

Given a discrete set of observation points of a deterministic 
or stochastic process i xCnO, n = 0,1..., N-1> we can fit them 
into the discussed AR model by Eqn (2.2.2) as : 

= E A ‘"i" **l * 

1 = 1 ^ 

= E A e-X”!" '’i "i" ' 

i=1 " 

P 

= E <5n,i (2.5.1) 

where A^^^ = Aj^ e^^i is the complex carrier amplitude of 
exponential . 

6 = ^1 
Tl y X 

Eqn (2.6.1) can be written as a matrix equation as follows 


— 

xC03 


6_ . 6_ ....6 


A 



0,1 0,2 o , p 


”c1 

xC13 

r=: 

^ ^ ^ <5 . « 6 


A 



1»1 1,2 ^1,p 


^c2 

• 


• • * 

I • 


• 


19 



xCN-13 


N-1,1 


N-1,2 


N-1,p cp 


NXp PX1 


or X = Z A 


(2.6.Z> 


Z IS a NXp rectangular matrix with N >> p. To solve eqn (2.6>2) we 
need to compute the inverse of Z. Me solve Eq (2.6.2) by using the 
pseudo inverse of Z. The singular value Decomposition (SVD) of 


matrix Z can be written as (173. 


X p = X N XpV*^ P X p (2.6.3) 

where Q has the left singular vectors of Z, 

D is a dia-gonal matrix containing the singular values of Z and 
V has the Right singular vectors of Z. Q and V are unitary 
matrices. 



where ad v. are column vectors and are the signular values 
of Z with 5: > 0. The psendoinverse of Z is given by 


Z* = V d' 

p X p 


_ - . fcl V* 


pxp pxN NxN 

The solution of Eq (2.6.2 is given as 

A = Z^ X 


(2.6.4) 


( 2 . 6 . 5 ) 



The magnitude of will give us the information about 
amplitude, and phase of will give the information about phase 
of the i^*' exponential carrier. 

2.7 CLOSE FREQUENCY CASE 

With establishment of the procedure of decomposing a signal 
as a sum complex exponential FM subsignals, we now proceed to a 
more complex case wherein the carrier frequencies are so close 
that the pattern of one overlaps with that of another. Clearly in 
such a situation, though the carrier and modulating frequencies 
can be ascertained from the PSD plots as explained before, the 

determination of modulation index, p , becomes more difficult as 
the magnitude of overlapping frequencies get altered from its true 
value. 

Under such ci rcumstances the advantage used in obtaining Eq 
(2.5.10) from Eq (2-5.9) no longer holds good. As such it is roust 
that pattern of each carrier and set of modulating frequencies on 
either side should be separated from one another, ie the affect 
of one waveform over another in DFT plot has to be removed. 

To overcome this problem, we bank on the properties of DFT, 
namely the real part of DFT is of even symmetry and imaginary part 
IS of odd symmetry C5D. So keeping the carrier frequency as centre 





we flip about other in forward direction or backward direction* 
thus creating the pattern about the carrier frequency. This 
pattern is then de4ucted from the overall signal thus removing its 
affect. Same process is carried on the residual signal till all 
patterns are separated. 

The process of flipping is started from one extreme carrier 
frequency as its frequency will be affected the least because of 
overlapping. Once the patterns are separated the value of the 
modulation index can be found by Eqn. (2.5.13). 

Due to dependence of resolution of DFT on the number of 
point of DFT taken, it is possible that at a moderate point DFT, 

two distinc-t peaks may appear as a single peak in the DFT 
magnitude plot. This type of possibility is predominantly high in 
side frequencies of two adjacent patterns. Under such 
circumstances using the above discussed flipping procedure will 
introduce a new frequency in the residual signal, which was not in 
the signal originally. 

To overcome this problem, we make an estimate of all possible 
frequencies whose combination can lead to the observed frequency, 
in the DFT magnitude plot. This is most easily done by joining the 
adjacent frequencies by linear interpolation. In short an envelope 
of the DFT magnitude plot for the frequencies of interest is 


obtained. 



The process of eliminating the affect of overlapping is then 
done by flipping as discussed above. After each flip, the spurious 
frequency appearing in the residual signal is removed and the 
process of flipping carried till all patterns are well resolved- 



CHAPTER - 3 


SIMULATION RESULTS 


3.7 Introducation 

To show the effectiveness of the model discussed in 
chapter 2, we do some simulation study with syntnesised data. A 
signal consisting of several compleK FM sub signals is sampled and 

the product function for different lags are computed. These are 

then fitted on an AR model whose coefficients are computed by Burg 
algorithm. As expected the PSD of the product function consists of 
sets of frequencies, with a high peak at twice the original 

carrier frequency with diminishing side peaks at true modulating 

frequency from the carrier. 

3.S SIGNAL DATA 1. 

The complex FM signal used for the purpose of simulation 
consists of two complex FM signals 

x(t)=E fl (n,l. + ^5^ sin rtjt) 

1=1 * 

where 

Al = = 800 n, = 30°, = .4, = 


70 n 



Ag = 1, Og = 1600 n, 10°, /?2 = 1, = 66 n 

Bandwidth <BW) of FM singal = 2{/^+1) fm 

ie BW. = 9S Hz 
1 

BWg = 132 Hz 

Nyquist Rate 1 = 1992 Hz 
Nyquist Rate 2 = 3728 Hz 

By sampling theorem sampling frequency (Fs) S: Nyquist Rate 
le Fs = 4000 Hz 

Satisfying the sampling theorem, the sampling frequency 
correponding to the period T=.25 m sec is chosen for uniform 
sampling grid. The samples taken are of the form 

>tCn3 =1 A e^ ^"i ^n + + /9. sin Njn),- n = 0,1,... 

i=1 ^ 

(3.2.1) 

The real and imaginary parts of the signal are plotted in Fig 
3.2.1 and Fig 3.2.2 respectively. 

3.2.1 SIMULATION RESULT 

Using the smaples generated by Eqn (3.2.1), the product 
function CK3 were calculated by Eqn (2.3.8) for K = (-100, -98, 
... -2, O, 2, .... 98, 100) a total of 101 points. The real and 
imaginary parts of P/:K3 are platted in Figs 3.2.3 and 3.2.4 
respectively. 


The PSD plot can then be obtained by utilizing the AR 



coefficients. In the PSD plot, if the earner and carrier plus 
modulating frequencies are not well resolved then we have to 
increase the model order to split them. However, if the AR model 
order chosen is too high then line splitting will occur ie- a 
single frequency may falsely split into two lines giving rise to 
spurious peaks. 


Figure 3.2.5 shows the effect of AR model order on the PSD 
plot. From the figure, for model order 50, the desired patterns 
are found centered at discrete frequencies 



2 ^2 = * 3 , 

2^3 “ 

Hence f^^ = .1, 

= - 15 . 

^3 = -2 

Clearly, the pattern 

.t = .15 = 

• 1 + .2 


term as discussed in Eqn (2.3.4), and as such of no interest to 
us. Hence the two carrier frequencies are f^^ = .1, f^^, = .2. 

Also from the same figure the side peaks give the modulating 
frequencies as 

^ m 1 = <^ ro 2 = •^'’2 

For evaluation of the modulation index le ft, we use the DFT 
method discussed in section 2.5. 

A 512 point DFT of the signal xCn3 is plotted in Fig 3.2.6. 
Calculating k (sample nos) corresponding to 1=0, 


± 1 , ±2 



using Eqn (S.5.7) and finding ratio Jo <^)/J1 </5), Jo (/?)/JS (/? ) 
etc using Eqn (E.5.13) for both patterns, in a tabular from we 
have , 


Set 

No 

l=-2 

K= 

1=-1 

K= 

1=0 

K= 

1 = +1 

K = 

L = +2 

K = 

Jo/JI Jo/J2 


1. 

— 

55.60 

51.2 

46.72 

- 

4.555 

.4 

2 . 

93.952 

98.176 

102.4 

106.624 

110.85 

1.751 6.219 

1 


The amplitudes and phase of the two signals are calculated 

using Eq (E.6.5) which give 

A. = 1, = 30.1° 

1 1 

Ag = 1, = 9.99° 

3.3 SIMULATION NO 2 

This simulation with the synthetic data deals with the 
procedure dealt in section 2.7. The complex FM signal used for 
the purpose consists of two complex FM signal with parameters as, 
A^ = 1, = 1400 n, = 30°, = .4, N^ = 70 n 

Ag ^ n, = 10°, = 1, Ng = 66 n 

Satisfying the sampling theorem, the sampling frequency 
corresponding to the period T=«25 m sec is chosen for uniform 
sampling grid. The samples taken are of the form of Eqn 3.2.1. 

Using the samples generated the product function P CK3 were 




2 


calculated by Eqn (2.2.&) by K = (—100, —98, .... 2, 0, 

98, 100). The product function <P„ Ck3) are then fitted in an AR 
model which uses Burg algorithm. From the PSD shown in figure 
3.3.1 the carrier and modulating frequencies are known. 


. = .175 
cl 

^2 

. = .014 
m1 

^m2 


A 1000 point DFT magnitude is plotted in Fig 3.3.2. By 
comparison with the patterns about the two carriers it is evident 
that the magnitude of frequencies between the two carriers have 
altered due to overlapping. 


Fig 3.3.3 shows the envelope obtained by linearly 
interpolating between the peaks. The interpolation is done on both 
the real and imaginary values of the complex DFT. Then they are 
mixed to find the magnitude. Working on real and imaginary part of 
DFT becomes necessary as the magnitude plot does not give any 
information about the manner of combination of two different 
complex values. 

First backward flip is carried about f^^ = .2. Pattern about 
this carrier frequency is well established by the right most curve 
in the Fig (3.3.4). This pattern is then deducted from the 
original signal envelope. The residual consists of the well 
desired pattern about frequency f^^ = .175, shown by the left 

curve of Fig (3.3.4). The spurious signal from the residual has 
been eliminated. 



Calculation of modulation indeK, from Fig (3.3.4) and Eqn 


CS.5-10) gives 




= 5.134 


which from tables give 


ft. 


.3a 


and for the other pattern 

which from tables gives 



<ft^^ 

~iftll 


6.61 


/?2 ^ .99 

The amplitudes and phase of the two signals can 
calculated using Eqn (2.6.5) 


be 




Fig. 3,2.1 Real Part of complex FM signal 



Fig. 3.2.2 Imaginary Part of complex FM signal 
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Fiq. 3.2.5a PSD s K = ± 100, Model order 50 
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CHAPTER -4 


SIHULATIOM RESULT WITH REAL TIME DATA 

4.1 INTRODUCTION 

As described in chapter 3 r the parameter estimation is found 
to be reasonably accurate for the model of sum of compl&j^ Fli 
signals. Also the procedure of estimation is shown to be 
performing well for the overlapping patterns case. The next step 
will be to try to fit some real time non-stationary data on this 
model. Speech signal is an example of a non— stationary signal 
showing frequency variations and hence was considered as a good 
candidate for this model. For this purpose, the model was fitted on 
one representative signal from each voiced and unvoiced speech 
signal -The representative signals were phonemes \k\ and \e\ 
respectively. 

4.S SIMULATION RESULTS WITH SPEECH DATAs 

, The speech data on which the model was fitted , consisted of 

1 

unvoiced \k\ and voiced \e\ .A total of 1500 sample points 
uniformly sampled at 16 KHz for each phonemes were available le 
the duration was slightly more than 90 ms. The Product Function was 
calculated for lag of — SOO to + 200 .These were fitted onto an AR 
Model and the PSD was computed . Starting with a low model order 

{ 

i of 10 and then increasing it, the correct model order was 


established. 



The PSD of unvoiced \k\ sound is shown in fig 4.2.1. Since 
the data are real valued and so are the Product Functions , the 
PSD plot is syiRinetr ical about ■F=0. The DFT magnitude plot is shown 
in Fig. 4.2.2 . From the PSD , taking aid of the DFT plot , the 
various FM carrier and modulating frequencies are determined . The 
spectrum consists of five FM sub— signals . From the DFT plot it is 
also clear that due to overlaping , many of the peaks have altered 
amplitudes. Also the peaks about the carriers are not at 

modulating frequency distance, indicating the poor resolutin of 
DFT. Dealing the DFT available , by method described in section 
2.7 , we find the envelope of the peaks as shown in Fig. 4. 2. 3 
Flipping about the estimated carrier frequency’s the pattern of 
each FM signal is separated. This is shwon in Fig. 4.2.4 .From 
each pattern, value of corresponding modulatin indeK is found. 
Using these frequencies and modulation indexes , we obtain the 
amplitudes and phases for the constituent FM signals. Using these 
paramenters we regenerate the consonant sound. Fig 4.2.5 shows the 
original and regenerated consonant sound. 

The PSD of voiced \e\ sound is shown in fig 4.2.6. Since the 
data are real valued and so are the Product Functions , the PSD 
plot IS symmetrical about ■P=0 . The DFT magnitude plot is shown in 
Fig. 4.2.7 . From the PSD , taking aid of the DFT plot , the 
various FM carrier and modulating frequencies are determined . The 
spectrum consists of two FM sub-signals . From the DFT plot it is 
also clear that due to overlaping , many of the peaks have altered 
amplitudes. Also the peaks about the carriers are at modulating 
frequency distance . Dealing the DFT available , by method 



described in section 2.7 , we find the envelope of the peaks as 
shown in F'i^.^.2.3. f"lippin9 about the estimated carrier 
frequency’s the pattern of each FM signal is separated. This is 
shwon in Fig. 4.2.9 . From each pattern, value of corresponding 
modulatin index is found. Using these frequencies and modulation 
indexes , we obtain the amplitudes and phases for the constituent 
FM signals. Using these parame-ters we regenerate the consonant 
sound. Fig 4.2.10 shows the original and regenerated consonant 
sound. 


The model has been shown to be suitable for modelling speech 
signals with fair degree of success. The model seems to perform 
better for the unvoiced speech case. This is expected as the 
modulation in amplitude is much less in unvoiced speech vis-a— vis 
voiced speech. In fact , by looking at the DFT plot, a fair idea 
can be obtained about the suitability of a particular signal to be 
modelled by this particular model. If the DFT plot consists of the 
desired patterns or approximates it, then the signal is a good 
candidate for this model. This is amply evident from the 
simulation cases dealt with. 
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CHAPTER-5 


CONCLUSION 

The model introduced in this thesis can be used to model 
signal data from a deterministic or a stochastic process , as a 
sum of several FM signals. With synthesized data, reasonably 
accurate estimation of parameters are obtained. 

The Fourier transform based techniques though limited in 
resolution, perform adequately under low SNR. AR modelling has the 
advantage of high resolution spectral estimation, though its 
performance is good at high SNR and deteriorate at low SNR .A 
combination of the two, using the capability of each makes the 
model parameter estimation feasible accurately. 

The model has been found particularly suitable for speech 
signals • Knowing the parameters of a signal the resulting 
parameters can be used as a feature vector for pattern 
recognition. The application of this model may include speech 
segmentation, speaker verification, speaker identification, speech 
recognition and language identification etc. 

However a much more detailed experimentation is required to 
fully exploit the capabilites of this model. Specifically a more 
general model consisting of a sum of AM— FM signals may be suitable 
to model any type of non— stationary signal. However, more work 



will be needed in this direction before the method can be 
standardised. 
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